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^ \ Abstract 
>. ■ 

Q>^ ■ The generalized Henon-Heiles system with an additional nonpolynomial term 

'nI" . is considered. In two nonintegrable cases new two-parameter solutions have been 

^ \ obtained in terms of elliptic functions. These solutions generalize the known one- 

\ parameter solutions. The singularity analysis shows that it is possible that three- 
^ ■ parameter single- valued solutions exist in these two nonintegrable cases. The knowl- 
^ ■ edge of the Laurent series solutions simplifies search of the elliptic solutions and 

p |. allows to automatize it. 

1 ■ 

a • 

S: 

> : 1 INTRODUCTION 

•i-H . 

;_i ■ Beginning from papers [1-3] , investigations of two-dimensional Hamiltonian systems with poly- 

nomial potentials attract large attention due to detect of the "dynamical chaos" phenomena. 
There is no method to find the multivalued general solution of a two-dimensional nonintegrable 
system in the analytic form. At the same time it is an actual problem to find single-valued 
special solutions in the analytic form, because the investigation of the solutions with some 
additional properties, for example, periodic solutions, plays an important role in the study of 
physical phenomena. Another problem is to pick out nonintegrable cases, in which single- valued 
special solutions can depend on maximal number of arbitrary parameters. 
The Henon-Heiles Hamiltonian [2]: 

H = \{x\ + vl + + 1/2) + a;2?/ - 1 

and its generalizations are one of the most actively studied two-dimensional Hamiltonians 
(see ^ and references therein). The generalized Henon-Heiles system is a model widely used 
in astronomy j3] and physics, for example, in gravitation j^l [7] . 
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One of lines of investigation of this system is the search for special solutions [8-13]. The 
general solutions in the analytic form are known only in the integrable cases [14-17], in other 
cases not only four-, but even three-parameter exact solutions have yet to be found. In ^2] 
new type of one-parameter elliptic solutions has been obtained. Such solutions exist only in 
integrable cases and two nonintegrable ones. In these nonintegrable cases there exist three- 
parameter Laurent-series solutions ^H], which generalize the Laurent series of one-parameter 
elliptic solutions. In this paper we find elliptic two-parameter solutions, which generalize solu- 
tions obtained in |12j . 

2 BASIC EQUATIONS 

The generalized Henon-Heiles system with an additional nonpolynomial term is described by 
the Hamiltonian 

1 c 

and the corresponding system of the motion equations: 

^tt = - Aix - 2xy + 
Vtt = - A2I/ - + Cy'^, 

where x^f- = ^ and y^^. = Ai, A2, and C are arbitrary numerical parameters. Note, that 
if A2 7^ 0, then one can put A2 = sign{X2) without the loss of generality. 

Due to the Painleve analysis [19-21] the following integrable cases have been found |22j : 

(i) C = -l, Ai = A2, 

(ii) C = —6, Ai, A2 arbitrary, 

(iii) C = -16, Ai = A2/I6. 

These integrable cases correspond precisely to the stationary flows of the only three inte- 
grable cases of the fifth-order polynomial nonlinear evolution equations of scale weight 7 (re- 
spectively the Sawada-Kotega, the fifth-order Korteweg-de Vries and the Kaup-Kupershmidt 
equations) [HI • 

In all above-mentioned cases system (2) is integrable at any value of /x. The function y, 
solution of system (2), satisfies the following fourth-order equation [TIH I12[ 1^: 

20(7 

yuu = (2C - 8)y,a - (4Ai + X2)yu + 2(0+ l)y^ + -^y' + (ACX, - 6A2)y' - 4Ai A2y - AH, (3) 

where H is the energy of the system. We note, that H is not an arbitrary parameter, but a 
function of initial data: ?/q, y^^, y^^^ and y^ftf The form of this function depends on /i: 

H--^r„2 ,,2N 3 /Ai \ 2 x„ „ ^ , i^2yot + ^Cy^ypt - yput)'^ + fi 
H - -(.0. + 1/0) - 3 .0 + + y.) iCy, - X.y, - y,,,) + ,^Cyi - X.y, - y,,) " 

This formula is correct only if Xq = Cy"^ — A2?/o 7^ 0- If ^0 ~ 0' what is possible only at 
/i = 0, then we can not express Xq^ through ?/q, ?/qj, t/q^^ and y^m, so H is not a function of the 
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initial data. If t/q^^^ = 2C?/q?/qj — \2y0t, then eq. (3) with an arbitrary H corresponds to system 
(2) with /i = 0, in opposite case eq. (3) does not correspond to system (2). 

To find a special solution of eq. (3) one can assume that y satisfies some more simple 
equation. For example, there exist solutions in terms of the Weierstrass elliptic functions, 
which satisfy the following equation: 

y^=Ay^ + By^ + Cy + V, (4) 

where A, B, C and V are some constants. 
The following generalization of eq. (4): 

= Ay' + By'/' + Cy' + Vy'" + Sy + G (5) 

gives new one-parameter solutions in two nonintegrable cases C = —16/5 and C = —4/3 
(Ai is an arbitrary number, A2 = 1). It is easy to show [E] that if S 7^ or P 7^ then Q = 0, 
therefore, substitution y = g' transforms eq. (4) into 

g''t=^[Ag'' + B0' + C0'' + Vg + £). (6) 

In [13] using the substitution y — ^ y — Pq a. new parameter Pq has been introduced and two- 
parameter solutions have been constructed for above-mentioned values of C and a few values 
of Ai (A2 = 1). Due to Painleve analysis local three-parameter solutions as the converging 
Laurent series have been found for an arbitrary Ai, A2 = 1 and /i = JHl- In the present paper 
we seek both the elliptic and the Laurent-series solutions for arbitrary values of Ai, A2 and fi. 



3 NEW SOLUTIONS 

Let us assume that solutions of eq. (3) in the neighborhood of singularity point to tend to 
infinity as y = cplt — to)^, where f3 and are some complex numbers. Of course, the real part 
of f3 has to be less then zero. From this assumption it follows [22] that f3 = —2. The Laurent 
series of solutions of eq. (6) begin with term proportional to (t — to)~^, so we seek solutions of 
eq. (3) as square polynomial: y = P2Q' + PiQ + Pq, where P2, Pi and Pq are arbitrary numbers, 
g is the general solution of eq. (6) with arbitrary coefficients A, B, C, V and £. Because of 
the function g = {g — -y) / is a solution of eq. (6) as well, we can put P2 = I and Pi = 
without loss the generality. 

Substituting y = g'^ + Pq in eq. (3), we obtain 

gtmg = - ^QtttQt - ^Qtt + 2(C - 4:)gug' + (2Po(C^ - 4) - 4Ai - X2)gttg + 
+ 2{3C - 2)q]q' + (2CPo - 4Ai - 8P0 - A2)f?? + ^C(f + 
+ {2C\i + lOCPo - 3A2)f?^ + 2(2AiCPo + WP^ - A1A2 - 3PoA2)f?' + 
+ ^CP^ + 2A1CP0' - 3P2A2 - 2A1A2P0 - 2P-. 
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The function p is a solution of eq. (6), hence, eq. (7) is equivalent to the following system: 



(3^ + 4) (-3^+ 2C) = 0, 
B{-21A + 9C - 16) = 0, 

96^CPo - 240AC - 192^Ai - 384^Po - 48^A2 - 

- 105i32 + 128CC - 192C + 128(:7Ai + 640CPo - 192A2 = 0, 
AOBCPo - 90AV - 65BC - SOBXi - WOBPq - 2OBX2 + 56CV - 6AV 
I6CCP0 - 36A£ - 21BV - 8& - 32CAi - 64CPo - 8A2C + 24C^ + 

_ _ , ^ 16S- 32A1A2 - 96P0A2 = 0, 

(5C + 8CP0 - I6A1 - 32Po - 4A2):P = 0, 
384i7 = -48C^ + 96CSPo + SSACX^P^ + 640Cp3 - 9p2 _ 

- 192^Ai - 384^Po - 48^A2 - 384A1A2P0 - 576A2Po^. 



0, 



(8) 



+ 64AiCPo + 160CP2 
lOBS 



System (8) has been solved by computer algebra software REDUCE 
If i3 7^ 0, then from two first equations of system (8) we obtain: 



C 



and A 



4 
3 



or 



C 



16 



and A = — 



32 
15' 



If = 0, then solutions with V ^ are also possible at C = —16 and C = —1, but only 
in integrable cases. The obtained solutions of eq. (3) depend on two parameters: energy H 
expressed through Pq and parameter to connected to homogeneity of time. 

Six solutions of system (8) correspond to each value of Pq. Two of them (with B = V = 0) 
generate solutions of eq. (4). Values of B and V, corresponding to other solutions, depend 
on Ai and A2 and are zero only at some relations between these parameters. We will consider 
only solutions with S 7^ or P 7^ 0. They are presented in Appendix. These solutions can 
be separated on pairs in such a way that solutions in one pair differ only in signs of B and V. 
Basic properties of the obtained solution are considered in this section. In the next section we 
analyze in detail solutions of system (8) for some values of Ai and A2. 

If the right-hand side of eq. (6) is a polynomial with multiple roots, then g and y can be 
expressed in terms of elementary functions. In opposite case y is an elliptic function |23] . 

It is simplicity itself that y(t) = g^it) + Pq = {—git))'^ + Pq, so, solutions of system (8) with 
opposite values of B and V generate identical solutions of eq. (3). From eq. (6) we obtain a 
polynomial equation for y{t): 



{yl - Aiy - Po)' - C{y - P,f - E{y - P,)f = (y - P,f{B{y - P,) + Vf. 



(9) 



The function g{t) can be expressed through the Weierstrass elliptic function p{t) Ch. 5]: 

ap{t -to) + b 



g{t - to 



cp{t - to) + d' 



(ad — be = 1), 



where to is an arbitrary parameter. Periods of p{t) and the constants a, b, c and d are determined 
by eq. (6). The function 

/ ap{t -to) + b 



yit - to) 



cp{t - to) + d 



(10) 
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is the fourth-order elhptic function. This function, as a solution of eq. (3), can have only the 
second-order poles, therefore, in the parallelogram of periods it has two poles with opposite 
residues. Solutions (10) differ from solutions of eq. (4), which are the second-order elliptic 
functions 

The function x{t) satisfies the first equation of system (2) with 



= ^C'Po' + (2\iC^ - y AsC) Po' + (2A^ - - 4A1A2C) + (2AiA^ - 

2 I /o\ \ C I C2 I /I \ Tj\ ID I nczj I \ c2 I Ti^c 



'111 



- 2X1CS - ACH + 3A2^ Pq + 2AiA2^ + + AX^H) Pq + 2£H + -\^£' + —V'E 

' ^ ' L Izo 



The trajectory of the motion can be derived from the second equation of system (2). Sub- 
stituting y^^, we obtain: 

x^ = {C- '^-A)y' + (3^Po - C - \)y - \{my + 3P - mP^)^y-P^- + 3^p2 - 2CPo). 

If B and T> take zero values we get simple algebraic trajectories. The full list of such 
trajectories is presented in jTT]. The parameter Pq is absent in these trajectory equations. 

One value of the energy H can correspond to no more than three values of Pq and, hence, no 
more than six different one-parameter solutions. Solutions (10) differ from solutions of eq. (4), 
which are the second-order elliptic functions |25j . 



4 A PARTICULAR CASE 



4.1 The form of solutions 



At C = —16/5, Ai = 1/9 and A2 = 1 one-parameter solutions (Pq = 0) have been considered 
in detail in our previous papers ^21 ^j. For these values of parameters solutions of system (8) 



are: 

~ 32 ~ -32 

1. ^ = --, i3 = 0, c = --Po-l, 

^ = 0, = - f Po^ - 2Po, ^ = I^^o+^^o 



4 ~ -17 
2. ^=-3, i3 = 0, C = -4Po--, 

~ , 34 20 2 ^ 17 , 2 230 

p = 0, E = -4P2 - — Po + E = - —Pi - —Pi + ^Po - 



33 3267 15 " 330 " 3267 323433 



32 R-+8^V15 ~_ 32 4 

15' ^-"Y^°"9' 



V ^ ±^Po, ~E ^ -'-^Pl - '-Po, H ^ 1^P3 - IPI, 
9 5 ° 9 15 ° 72 °' 
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5-6. i=-|, l^.±^^(26928P„ + 8125), 



"8415 ' 5 ^ 1683 7553304' 

~ 32 „ 1748 IQ. 7291 6426875 „ 17551324375 



5 1683 15 ^ 13464 ^ 181279296 9762977765376 

If the right-hand side of eq. (6) is a polynomial with multiple roots, then the function y can 
be expressed in terms of elementary functions. For example, at Pq = substitution of solutions 
3-4 into eq. (5) gives 

y = -— — — r-T' (12) 



3(l-3sin(^)) 



where to is an arbitrary constant. 

From (11) we obtain the following values of /x: 

1. fi = Q, 

9 /, = 160 p3 I 680 p2 _ 800 p _ 7000 

1089 "T 11979 1185921 1056655611' 

O 'i. — g-TQ -|- g^J-Q -|- 729-* 0' 

e. _ f. 52 p4 81640 p3 4458460825 p2 539878421875 p 728473377734375 

^ ^- 561 944163 152546527584 128367902961936 6703885364284145664' 

4.2 Motion trajectories 

Let us consider the equations of the motion trajectories at C = —16/5 and A = 1/9. In the 
case of the solutions with B = V = the trajectory equation can be reduced either to = 
(solution 1), or to 

. 6 / 20\2 50 

X + - [y + — ] = . 13 

5V 99/ 1089 ^ ^ 

In the last case (solution 2) the motion trajectory is an ellipse. Note, however, that the real 
motion does not necessarily affect the whole ellipse: it depends on two arbitrary parameters. 
The energy H can be considered as one of them. 

In the case of solutions 3-4 the trajectory equation is the following: 

+ ly)' + §fiy - Po)i^y + Po)' = o. (m) 

If Pq = (see (12)), the equation for one of the trajectory branches entirely coincides with 
the equation obtained in The condition y < is always required for the existence of real 
motion along these trajectories. Formula (9) describes precisely such a solution. For solutions 
5-6 the trajectory equation has the same form as for solutions 3-4. 



5 THREE-PARAMETER SOLUTIONS 



The Ablowitz-Ramani-Segur algorithm of the Painleve test ^20j is very useful for obtaining the 
solutions as formal Laurent series. Let the behavior of a solution in the neighborhood of the 
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singularity point to be algebraic, i.e., x and y tend to infinity as some powers: x = a^it — to)"- 
and y = hpit — to)^, where a, /3, and hp are some constants. If a and (3 are negative 
integer numbers, then substituting the Laurent series expansions one can transform nonlinear 
differential equations into a system of linear algebraic equations on coefficients of Laurent 
series. If a single-valued solution depends on more than one arbitrary parameters then some 
coefficients of its Laurent series have to be arbitrary and the corresponding systems have to 
have zero determinants. The numbers of such systems (named resonances or Kovalevskaya 
exponents) can be determined due to the Painleve test. 

Two possible dominant behaviors and resonance structures of solutions of the generalized 
Henon-Heiles system (22 12ZI and eq. (3) are presented in the Table. 



Case 1 


Case 2: (3 < 3?e(a) 


a = -2, 






1 ± 7l - 48/C 

a = , 

2 




/5 = -2, 






/? = -2, 




a„ = ±3^2 4 


-c, 




o-a = Ci (arbitrary). 




h = -3, 






bp- 




r = -l, 6, - 


^ yi - 24(1 + 
2 


C) 


r = -l, 0, 6, Tsjl- 


48 


r, = -1, 10, 


5^ yi-24(l 
2 2 


+ C) 


r4 = -1, 5, 5- 


48 ^ 1 48 
C'^ + f -C 



The values of r denote resonances: r = — 1 corresponds to arbitrary parameter to; r = 
(in the Case 2) corresponds to arbitrary parameter Ci- Other values of r determine powers 
of t, to be exact, f^^"' for x and t^"*''" for y, at which new arbitrary parameters can appear as 
solutions of the linear systems with zero determinant. Note, that the dominant behaviour and 
the resonance structure depend only on C. 

It is necessary for the integrability of system (2) that all values of r be integer and that all 
systems with zero determinants have solutions for any values of the free parameters entering 
these systems. This is possible only in the integrable cases (i)-(iii). 

For the search for special solutions, it is interesting to consider such values of C, for which 
r are integer numbers either only in Case 1 or only in Case 2. If there exist a negative integer 
resonance, different from r = —1, then such Laurent series expansion corresponds rather to 
special than general solution [221 ■ We demand that all values of r, but one, are nonnegative 
integer numbers and all these values are different. From these conditions we obtain the following 
values of C: C = -1 and C = -4/3 {Case 1), or C = -16/5, C = -6 and C = -16 {Case 

/ 1—48 /C 

2, a = — ^^"2 ), and also C = —2, in which these two Cases coincide. It is remarkable that 

only for these values of C there exist solutions of system (8) with i3 7^ or P 7^ 0. 

Let us consider the possibility of existence of the single- valued three-parameter solutions in 
all these cases. To obtain the result for an arbitrary value of ^, we consider eq. (3) with an 
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arbitrary H. Note, that the values of resonances obtained from eq. (3) (in the Table they are 
signified as r^) are different from r, but we obtain the same result: condition that all values of 
r4, but r4 = —1, are nonnegative integer numbers gives the same values of C. 

At C = —2 we have a contradiction: r4 = 0, but 6_2 is not arbitrary parameter: 6_2 = —3. 
This is the consequence of the fact that, contrary to our assumption, the behaviour of the 
general solution in the neighborhood of a singular point is not algebraic, because its dominant 
term includes logarithm [22] • At C = —6 and any value of other parameters the exact four- 
parameter solutions are known. In cases C = —1 and C = —16 the substitution of an unknown 
function as the Laurent series leads to the conditions Ai = A2 or Ai = A2/I6 accordingly. Hence, 
in nonintegrable cases three-parameter local solutions have to include logarithmic terms. Single- 
valued three-parameter solutions can exist only in two above-mentioned nonintegrable cases: 
C = -16/5 and C = -4/3. 

Using the method of construction of the Laurent series solutions for nonlinear differential 
equations describing in [TH], we obtain single-valued local solutions of eq. (3) both at C = —16/5 
and at C = —4/3. Values of other parameters are arbitrary. 

At C = —4/3 these solutions are: 



1 , 1 29,0 1, 



4 

( ^x^b' 1 + —Xl + -A^ + 1 + 62 - -A26' 1 - — A1A2 ) b' + 

V12 90 ^ 40 ^ 2592 "i ^ ^ 16 60 ' 



A2+(-6^ + -A, 



4 
131 



A2 b-it + b2f 



(15) 



There exist four possible values of the parameter 6_ 



6_i = ±^ 



or 



6_, =± 



105A2 


- 140Ai + . 


y7(1216A? - 


I824A1A2 + 783Ai) 


i 385 


IO5A2 


- 140Ai - ^ 


/7(1216Af - 


I824A1A2 + 783Ai) 


385 



The parameters 62 and bs, coefficients at and correspondingly, are arbitrary. The energy 
H enters in coefficients beginning from 64. 

At C = —16/5 we obtain the following solutions: 



15 r 5 , 62-2 / 5 , 632-2 25 , , , 

y = ^ + 6-1 A2 H 6 1 + — Ai H r 1 A2 b_it 

^ 8t"2 32 45 V12 225 192 ' 

/29 ~2 1 ,2 29 , -2 102272-4 \ 2 r s 

+ —Xib^ 1 A2 X2b 1 + 1 U + + • • • 

Vl5 128 ^ 48 10125 "V ^ ^ ^ 



(16) 



with 



b-i 



±3 
41888 



6872250A2 - 21991200Ai + 52360 J71680Af - 448OOA1A2 + 13545Ai 



or 



±3 
41888 



6872250A2 - 21991200Ai - 52360V71680A? - 448OOA1A2 + 13545Ai. 
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The coefficients 63 and are arbitrary parameters. Beginning from 64 some coefficients 
include the energy H. So, the obtained local solutions depend on four independent parameters: 
to, H and two coefficients (62 and or 63 and fog)- 

We have found local single- valued solutions. Of course, existence of local single- valued solu- 
tions is necessary, but not sufficient condition to exist global ones, because solutions, which are 
single- valued in the neighborhood of one singularity point, can be multivalued in the neighbor- 
hood of another singularity point. So, we can only assume that global three-parameter solutions 
are single-valued. If we assume this and moreover that these solutions are elliptic functions 
(or some degenerations of them), then we can seek them as solutions of some polynomial ffist 
order equations. There are a few methods to construct such solutions [SI CHI I2H1 12^] • Using 
these methods one represents a solution of a nonlinear ordinary differential equation (ODE) as 
the finite Taylor or Laurent series of elliptic functions or degenerate elliptic functions, for ex- 
ample, tanh(t). Similar method is applied in this paper to find two-parameter solutions. These 
methods use results of the Painleve test, but don't use the obtained Laurent-series solutions. 
In 2003 R. Conte and M. Musette jHO] have proposed the method, which uses such solutions. 

The classical theorem, which was established by Briot and Bouquet |3^, proves that if the 
general solution of a polynomial autonomous ffist order ODE is single- valued, then this solution 
is either an elliptic function, or a rational function of e"'^, 7 being some constant, or a rational 
function of x. Note that the third case is a degeneracy of the second one, which in its turn is a 
degeneracy of the ffist one. It has been proved by Painleve JH] that the necessary form of the 
polynomial autonomous ffist order ODE with the single-valued general solution is 

m 2m~2k 

E E h^ky'yt = o, K^ = i, (17) 

A;=0 j=0 

in which m is a positive integer number and hj^ are constants. 

Rather than to substitute eq. (17) in some nonintegrable system, one can substitute the 
Laurent series of unknown special solutions, for example, (15) or (16) in eq. (17) and obtain 
a system, which is linear in hjj^ and nonlinear in the parameters, including in the Laurent 
coefficients |30, . There are a few computer algebra algorithms which allow to obtain this system 
from the given Laurent series. Moreover it is possible to exclude all hjj^ from this system and 
obtain a nonlinear system in parameters of nonintegrable system and free parameters from the 
Laurent series. The main preference of this method is that the number of unknowns in the 
resulting nonlinear algebraic system does not depend on number of coefficients of the ffist order 
equation. For example, eq. (17) with m = 8 includes 60 unknowns hj^^, and it is not possible 
use the traditional way to find similar solutions. Using this method we always obtain nonlinear 
system in 5 variables: Ai, A2, H and two arbitrary coefficients of the Laurent- series solutions. 
We hope that this method allows us to find three-parameter global solutions. 

6 Conclusions 

Two nonintegrable cases (C = —16/5 or C = —4/3, Ai, A2 and /x are arbitrary) of the general- 
ized Henon-Heiles system with the nonpolynomial term have been considered. To avoid prob- 
lems with the nonpolynomial term we have transformed system into the fourth-order equation. 
Two-parameter elliptic solutions for this equation have been found in both above-mentioned 



9 



cases. Two different solutions correspond to eacli pair of parameter values. The Painleve test 
does not show any obstacle to the existence of three-parameter single-valued solutions, so, the 
probability to find exact, for example elliptic, three-parameter solutions, that generalize the 
obtained solutions, is high. 

S.Yu.V. is grateful to F. Calogero, R. Conte, V. F. Edneral and A. K. Pogrebkov for 
valuable discussions. This work has been supported by Russian Federation President's Grants 
NSh-1685.2003.2 and NSh-1450.2003.2 and by the grant of the scientific Program "Universities 
of Russia" . 



APPENDIX 

In two nonintegrable cases (C = —16/5 and C — —4/3) for arbitrary Ai and A2 we obtain that 
six solutions of system (8) correspond to each value of Pq. Two of them (with B = V — 
generate solutions of eq. (4). Other solutions of system (8) can be separated on pairs such as 
each pair of solutions corresponds to one two-parameter function y = + Pq, where g satisfies 
eq. (6) with the following values of coefficients: 

5' 
15' 



^ v^TT22(1120Ai + 4I888P0 + 655g + 6I95A2) ^^1(^1, A2, Pq 



29373960(3600Af - II2OA1P0 - 2425AiA2 - 20944Po^ - 6I95A2P0 + 225Ai) ' 

2^°A 32 4 112 

187^^" 5 °+ 1309'^^~ 187^^' 



88320 2 480, „ 885 , ^ 153375. 



244783 ' 187 244783 ^ 244783 

32p2, 8 224 685 168855 2 

" + 1309^°'^^ " 187'^^^'^ " 3916528'^^'^^ + 3916528-^^' 

_ 11516270 3 8740 2p _ 3296515 50336425 2 , 

^ " "45774421'^^ + 34969 '^^ ^ ~ 2563367576'^^'^'^ + 183097684 -^^-^^ + 

258 2 8209 76915 12202395 

+ 187^^^° ~ 1958264^^^°'^'' + 279752^^-^^^° + 82027762432^^^^'^'' " 
131879855 , ,0 43 103 , ^2 8881 



_ Al A9 — PnSn ~\~ AoPn H XoPnSa — 

11718251776 ^ 13090 ° ^ 1496 ° 31332224 * 
71205 2p _ 12990165 ^ _ 168661575 3 16^3 
4476032 ^ ° ~ 1312444198912 ^ ^ ~ 187492028416 ^ ^ 15 ° ' 
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4-4 
(7 = -- A = -- 
^ 3' ^ 3' 



V330(952Ai - 616Po + 13i?, - 945A2)^F2(Ai, As, Pq) 
^ ~ 38115(432A? + 952AiPo - 29IA1A2 - 308Po^ - 945PoA2 + 27X 

C = Ai — 4Pn — — Rn — — Ao, 

33 " 66 ^ 22 ' 



2\ ' 



~_ 3394 2,54 1123 5897 

17 2 31 349 1223 13005 2 

3 308 ° 44 ^ ° 27104 ^ 3872 2' 

552922 29801 2 „ 173605,2^ 778033, 2 ^ ^ 185, . 

H = Ar AfPo H AfP„ H Af A2A2 AiPn + 

83853 ^ 2541 ^ " 2347884 ^ ^ 74536 ^ 66 ° 

20328 1 9 ^ 6776 ° 12522048 ^ 596288 ^ 1232 ° 

5 3 865 _ ^2 3065 _ ^„ 225909, 2d 2733 _2d 57699 ,3 

~l~ ~Pn ~l" — Z — AoPn — AoPn-flo — AnPn -|- z^Ao-Tlf, -|- — Ao, 

2 ° 176 ° 54208 ^ 54208 ^ 260876 ^ ^ 74536 ^' 

where 

Pi(Ai, A2, Po) = 39474176000A? + 122782105600A?Po - 104358400A^5, - 

- 17822336000A?A2 + 210552545280AiPo2 - 680261120AiPoSg - 10941145600AiA2Po - 

- 41066800AiA2Sg + 8305290000AiA2 - 501315584PoX - 65797670400A2Po + 



+ 55920480Po^g + 1611640800A2Po + 2884725A2^5 - 468507375A 



3 

Sg = ±^35(2048Af - I28OA1A2 + 387Ai), 

P2(Ai, A2, Po) = 2099776A? - 497728A?Po - 20008A?Pg - 4911144A?A2 + 948640AiPo^ + 
+ 19096AiPoPg + I458O72A1A2P0 + 37173AiA2Pg + 3943233AiA^ + Q77QP^Rq - 
- 7II48OA2P0 - 9240A2PoPg - 615384A^Po - 13581A^Pg - 1006425A^, 



Rq = ±A/7(1216Af - I824A1A2 + 783A 
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